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ABSTRACT 

We present a technique that combines Zeeman Doppler imaging (ZDI) principles 
with a potential field mapping prescription in order to gain more information about the 
surface field topology of rapid rotators. This technique is an improvement on standard 
ZDI, which can sometimes suffer from the suppression of one vector component due 
to the effects of stellar inclination, poor phase coverage or lack of flux from dark areas 
on the surface. Defining a relationship beween the different vector components allows 
information from one component to compensate for reduced information in another. 
We present simulations demonstrating the capability of this technique and discuss its 
prospects. 
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1 INTRODUCTION 



et 



al. 2000|; [Petit et al. 2000| ; [Ponati et al. 2000| ). As |CoL 



The so lar-st e llai" analogy is commonly invoked to explain 
stellar activity phenomena in terms of solar features such 
as prominences and flares (Radick 1991). A solar-type dy- 
namo mechanism is also thought to operate in other late- 
type stars. However it is unclear if this exact same mecha- 
nism will operate in stars covering a wide range of convec- 
tion zotie dSpthfi, l-otati6ti ta.t6s atid Marisfis. Thfotigh thd 
Mt Wil 



licr Cameron et al. (2000) report, these measurements lend 
support to dilTerential rotation models in which rapid rota- 
tors are found to have a solar-typ e differential rotation pat- 
tern (Kitchatinov & Riidiger 1999). Doppler maps of rapidly 



rotating stars show flux patterns which are very different to 
those seen on the Sun, with polar and/or high- latitude struc- 
ture often co-existing with low- latitude flux (De Luca, Fan 



II&.K 



[TcTr iOLd- 



suivy^, wliiLli lias inuniLuied sLblld 
tion periods and activity cycles using chromospheric emis- 
sion variations spanning a period of over 30 years, we know 
that many cool stars also display regular solar-type activ- 
ity cycles. However, other stars in the survey show irregular 



Saar 1997 



Strassmeier 1996 ). This is hard to reconcile 



with the solar distribution where sunspots and therefore the 
greatest concentrations of flux tend to be confined between 
±30° latitude. Does this indicate diff erent dynamo modes 



are excited in these stars ? Models by Granzer et al. (2000) 



variations and some show no cyclic patl 


erns over the time 


span of the observations (iDonahue 1996 


; Baliunas & Soon 


1995). Saar & Brandenburg (1999)1 and 


Brandenburg, Saar 


& Turpin (1998) find that the type of dynamo activity ap- 



and schiissler et al. (1996) explain the presence of mid-to- 
high latitude structure on young, rapidly rotating late-type 
stars in terms of the combined effects of increased Coriolis 
forces and deeper convection zones working to "pull" the 



pears to change with rotation rate and age in a study using 
data from the Mt Wilson survey as well as other photometric 
studies of cool stars. 

Doppler imaging techniques can map Teft fiux distri- 
butions on the surfaces of rapidly rotating cool stars JCof 
lier Cameron 1992: Piskunov, Tuominen fc Vilhu 1990; Rice 



flux closer to the poles. Even though [Cranzer et al. (2000) 



Wehlau fc Khokhlova 1989| ; [Vogt 198^ ). Starspot maps 
of cool stars taken a few rotation cycles apart allow us 
to determine the differential rotation rate on rapid rota- 



tors accurately (Ponati fc Collier Cameron 199 



Barnes 



flnd that both equatorial and polar fiux can exist on T Tauri 
stars that have very deep convection zones, in terms of flux 
"slipping" over the pole, they find the emergence of low- 
latitude flux on K dwarf surfaces difficult to explain. 

The technique of Zeeman Doppler imaging (ZDI) allows 
us to map the magnetic field distributions on the surfaces 
of rapid rotators usi ng high resolution circularly polarized 
spectra ( Semel 1989| ). Magnetic field maps of the subgiant 
component of the RS CVn binary, HR1099 {Prot = 2.8d), 
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and the KO dwarfs, AB Dor {Prot = 0.5d), L Q Hya (P, 



0.00 



1.6d) have been presented in several papers (Donati fc Col- 



her Cameron 1997 



Donati et al. 1999 



Donati 199£). These 



maps all show patterns for which there is no solar coun- 
terpart: strong radial and azimuthal flux covers all visible 
latitudes. On AB Dor, strong unidirectional azimuthal field 
encircling the pole is consistently recovered over a period of 
three years. This is very different to what is seen on the Sun, 
where hardly any azimuthal flux is observed at the surface 
and mean radial flelds are much smaller than those observed 
in ZD maps. 

While ZDI is an important tool in measuring the sur- 
face flux on stars, some questions about the authenticity 
of the features in these Zeeman Doppler maps remain. ZDI 
makes no assumptions about the field at the stellar surface; 
the radial, azimuthal and meridional ve c tors are mapped 



1997 |; |Ifrown et al. 1991[ ). For the true stellar field, how- 
ever, the physics of the stellar interior and surface will de- 
termine the relationship between the field components. Ad- 
ditionally, poor phase coverage is found to lead to an in- 
creased amount of cross-talk between radial, meridional and 



com pletely independently ( Hussain 199!: ; Donati & Brown 



azi muth al field components in ZDI maps (Donati & Brown 
1997j)7J) ue to the pre sence of starspots, t hese maps are also 



fiux-ce nsored. Indeed Donati et al. (1999) and Donati fc Col ^ 
lier Ca |neron (1997) find that the low surface brightness 



regions on AB Dor appear to correlate with the strongest 
regions of radial magnetic flux. For these reasons Zeeman 
Doppler (ZD) maps may not present an accurate picture of 
the surface held on these stars. 

By mapping potential fields on the surfaces of stars, we 
can evaluate more physically realistic models of the surface 
fiux distribution. While we do not necessarily expect the 
field to be potential loc ally, this assumption s hould be ade- 



quate on a global scale (Demoulin et al. 1993). Jardine et al _ 
(1999)|Jevaluate the consistency of Zeeman Doppler maps 
for AB Dor with a potential field using maps obtained over 
three years. Although they find that the potential held can- 
not reproduce features such as the unidirectional azimuthal 
fiux found at high latitudes, they do find a good correlation 
between the ZD azimuthal map and potential field conflgu- 
rations produced by extrapolating from the ZD radial map. 

The technique presented in this paper reconstructs po- 
tential field distributions directly from observed circularly 

pnlari'zpH profilpt; Thrt; t.prhniqiip allnwt; nt; tn prnHnrp Hif- 



ferent donfiaurations which are more phvsicallv realistic and 
which provide more information about the surface topology. 



2 TE E TECHNIQUE 




0.4 0.6 
Filling factor 

Figure 1. The power in each image pixel, i, (in each Im mode) is 
allowed to vary between —Bmax (i.e. filling factor, fi = 0.0) and 
+Bmax {fi = 1.0), where Bmai = lOOOG. The entropy is at a 
maximum value for a zero magnetic field {fi = 0.5). In the absence 
of information the code will push to the maximum possible value 
of entropy. Hence it weights equally for both positive and negative 
polarities. 
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Here Pim{d) is the associated Legendre function at each 
latitude, 0. 

We map the components, aim and bim - They are allowed 
to vary until they produce surface field vector distributions, 
Br, B^ and Be, that fit the observed polarized profiles. The 
procedure by which the circularly polarized dataset is pre- 
dict ed u sing a Br, B^ and Be image is described in Sec- 
tion 2.3. As with conventional Doppler imaging, this prob- 
lem is ill-posed and many different aim and bim image conflg- 
urations can fit an observed dataset to the required level of 
. A regularizing function is thus required to reach a unique 



solution and we use maximum entropy in this capacity ( Gull 



Skilling 1984 ; Skilling fc Bryan 1984). The maximum en- 



tropy method has been used in several Doppler imaging 
codes and enables us to produce a unique image which has 
the minimum amount of information that is required to fit 
the observed dataset within a certain level of misfit (as mea- 
sured by the statistic) ( golIi£r_CameronjJ 



ruh 1992; Piskunov, Tuominen & Vilhu 199C; Rice, Wehlau 



2.1 Modelling potential fields 

The surface magnetic field, B, is defined such that B=— V^*. 
The condition for a potential field, VxB=0, is thus satisfied 
and V.B=0 reduces to Laplace's equation: 



Khokhlova 198^ ) 



0. 



(1) 



At the surface, the radial field, Br, azimuthal field, B^, 
and meridional fleld, Be, can be expressed in terms of spher- 
ical harmonics: 



2.2 Maximum entropy 

We store magnetic field strength in terms of a filling factor 
model. The power in each mode is defined as a fraction of 
a maximum value of 1000 G. This is taken to be the upper 
limit for which the weak field approximation can still be 
applied (see Section 2.3). Magnetic fiux values between the 
lower and upper limits, -1000 G and 1000 G, are mapped 
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Figure 2. A schematic diagram illustrating how circularly polarized spectral signatures from magnetic field spots vary depending on the 
field vector. The contributions from radial, azimuthal and meridional spots to the summed Stokes V spectral signatures are represented 
by different line types (as shown in the key). 



onto flux filling factors, fi, between 0.0 to 1.0 (hence fi — 
0.5 = 0.0 G). 

For each mode, i, the entropy, S{f), is then defined in 
terms of fiux filling factor, fi, and image weights, Wi: 



S{f) = -e™.[/., 



4+(i-i^)iogi34] 



(5) 



The default image value, d = 0.5 (i.e. G). Fig. |l| shows 
how this form of entropy pulls equally in both directions to 
zero field images. 

On assigning equal weights to all Im modes, we find 
that the reconstructed images are preferentially pushed to 
the highest, most complex distributions. However we want 
to obtain the simplest possible potential field image which 
fits the observed spectral dataset. In order to obtain this 
reconstruction, we assign weights according to the number 
of nodes at the surface for each mode (i.e. each Im combi- 
nation) : 



W{l,m] 



= [(2lm| + l)(Z-|m| + l)]^ 



(6) 



The choice of n determines how strongly the image is 
pulled to the minimum number of nodes: n = 1 means that 
the weights are proportional to the number of nodes on the 
surface; n = 2 represents the surface area density of nodes. 

It should be noted that this weighting scheme simplifies 
the amount of information we reconstruct about the stellar 
surface. The choice of n in Eqn. ^ does not affect the sim- 
ple reconstructions but is found to produce more accurate 
reconstructions in the case of the more complex topologies. 
This is discussed in more detail in the next section. Using the 
solar-stellar analogy, it is likely that the actual distribution 
of magnetic regions is much more complex. The weighting 



scheme used here allows us to test the code and obtain a 
better idea about fundamental limitations inherent in the 
scheme. When dealing with real stellar spectra, a more re- 
alistic weighting scheme may be to push the reconstructed 
image towards a power-law distribution. The detailed effect 
of different weighting schemes on real stellar spectra will be 
the subject of a later paper. 

2.3 Modelling circular polarization 

The procedure used to calculat e predicted polarized line pro- 

Briefiy, 
This in- 



files is described more full y in Hussain ct al. (2000) 



Donati fc Brown (1997) 



the model used is that of 
volves making the following assumptions about the intrinsic 
line profile: 

• that it is constant over the surface; 

• it can be modelled using a Gaussian; 

• the weak field approximation is valid. 

Within the weak field regime, when Zeeman broadening 
is not the dominant line broadening mechanism, the Stokes 
V profile contribution, Vi, at a, wavelength. A, from a mag- 
netic region, s, is found to behave as follows: 



Vs{v) oc g\ 



dlsjv) 
dv 



(7) 



Here g is the average effective Lande factor over all the con- 
sidered transitions; Is is the unpolarized intensity line pro- 
file contribution from region, s; and i; is the velocity shift 
from the rest wavelength. The weak field approximation is 
found to be valid for mean (least squares de convolved) pro- 
files up to magnetic field strengths of 1 kG ( Bray & Lough- 
head 196E| ). 
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Through rotational broadening, flux contributions from 
different regions on the stellar surface are separated in veloc- 
ity space. As with conventional Doppler imaging, the veloc- 
ity excursions of the line profiles indicate the region of the 
stellar surface where the signature originates. High-latitude 
field signatures are confined to the centre of the line pro- 
file and distortions at the line profile wings are caused by 
magnetic fields in the low latitude regions. Circularly polar- 
ized (Stokes V) profiles are only sensitive to the line-of-sight 
component of the magnetic field. Fig. ^ is a schematic dia- 
gram showing how magnetic field spots with similar parame- 
ters but different vector orientations have differing Stokes V 
contributions depending on their position on the projected 
stellar disk. The Stokes V signature from an azimuthal mag- 
netic field spot switches sign as the spot crosses the centre 
of the projected stellar disk. Stokes V signatures from radial 
field spots do not show any switch, just scaling in amplitude 
with limb darkening. Similarly, low-latitude meridional field 
spot signatures do not switch sign. However, as the line-of- 
sight projection of low-latitude meridional vectors is sup- 
pressed, the circularly polarized signature is also reduced. 
Hence amplitude variations in a time-series of circularly po- 
larized spectra indicate both the strength of the magnetic 
flux at the surface as well as the vector orientation of the 
field. 

For an exposure at phase fe, the Stokes V specific in- 
tensity contribution, Vs, for a surface element pixel, s, to a 
data point in velocity bin, j, is determined as shown below 
for each fleld orientation: 



B„ 



-[V{A\sjk,fJ-sk)cOS Osk]. 



(8) 



Here Ba is the magnetic field value in pixel, s, at a particular 
field orientation and Bmax = 1000 G. V{AXsjk) represents 
the Stokes V specific intensity contribution for a 1000 G field 
pixel that has been Doppler-shifted by the instantaneous 
line-of-sight velocity of the pixel. The term, cos 9sk, is the 
projection into the line-of-sight of image pixel, s, at phase, 
k. 



3 SIMULATIONS AND RECONSTRUCTIONS 



We use a similar approach to Donati & Brown (1997) 



and 

demonstrate the capabilities of this technique by compar- 
ing potential field and Zeeman Doppler reconstructions for 
several different input datasets. 

All reconstructions presented here have been fitted to a 
reduced = l-O- For both sets of images, we used stellar 
parameters based on the active KO dwarf, AB Dor: P,ot = 
0.51479 d; inclination angle, i — 60°; Ve sin i— 90kms~^; 
'^rad^ Okms-\ 



3.1 Simple images 

The line and image parameters used for this test reconstruc- 
tion are listed below. 

Image: ai,i, = 300.0 G; = -300.0 G; 63,2 = 250 G 
Line: Ao = 5000 A, S/N = l.OE-l-05, spectral resolu- 
tion — 3 kms~^ 



Table 1. Reduced chi-squared values comparing input and re- 
constructed images. The first column lists which images they 
correspond to, the remaining three columns indicate the level of 
agreement between reconstructed images for each field orienta- 
tion. Clearly the potential field reconstructions show better agree- 
ment with the input images compared to the Zeeman Doppler 
reconstructions. 



Image 








Xmer 


Full Inp - 


ZD 


40.6 


57.2 


15.5 


Pull Inp - 


PF 


36.1 


12.4 


6.6 


50% Inp - 


ZD 


50.3 


54.2 


17.2 


50% Inp - 


PP 


39.5 


13.8 


9.3 



The Z D maps were p roduced using the ZDI code de- 



scribed in fHussain (1999) 



The input surface field, recon- 
structed potential field and Zeeman Doppler maps are shown 
in Fig. 1^. The model Stokes V spectra and fits from the po- 
tential field reconstructions are plotted in Fig. ^. The field 
in the visible hemisphere (—30° < 6 < 90°) is considerably 
smeared out in both the potential field and ZD reconstruc- 
tions. There is also clearly some cross-talk between the high- 
latitude meridional and azimuthal components in the ZD 
maps. This is because circularly polarized signatures from 
both high-latitude azimuthal and meridional fields change 
sign over a full rotation phase. This leads to ambiguity in the 
reconstructions. Low-latitude radial and meridional fields 
suffer from a similar ambiguity as the line-of-sight vector 
component from these regions does not change sign thus 
leading to cross-talk. As the potential field maps show, the 
additional constraint of a potential field leads to less am- 
biguity in the reconstructions. While the polarities are re- 
constructed correctly for the potential meridional field map, 
the flux strengths are greatly suppressed. This is because 
the meridional fi eld vector is a l most completely suppr essed 
at low latitudes (Hussain 1999; Donati & Brown 1997) and 
hence the information needed to reconstruct them is largely 
absent in the circularly polarized proflles 

The agreement between input and reconstructed images 
has been evaluated by calculating the reduced value be- 
tween each input and reconstructed map. The flux values 
for each magnetic fleld map were first normalised according 
to pixel area, such that the sum of all pixel areas is equal 
to 47r. Where A''= total number of pixels in each map, Ji = 
corresponding pixel in the input map and i?i=reconstructed 
map pixel; the reduced value is then evaluated is deflned 
as X^/i^ ^ 1) where: 



JV 



R^) 



As shown in Table |i| potential field reconstructions in 
all three field vectors show a better level of agreement with 
the input maps compared to the ZD maps. The reduction 
in cross-talk between surface field vectors is refiected by the 
much lower reduced values in the PF meridional and 
azimuthal field maps. 
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Input radial Image 



PF radial field map 



ZD radial field map 



fa, : I 2^ 
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longitude 
ZD azimutiioi field mop 
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longitude 
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50 180 210 240 270 300 330 360 

longitude 



Figure 3. Input surface radial, azimutlial and meridional field maps calculated using ai_i = 300.0 G; f)i,i = —300.0 G; 63,2 = 250 G. 
Grayscales represent the magnetic flux strength and the tick marks denote the observation phases. Potential surface radial, azimuthal 
and meridional field maps are plotted in the middle column. Zeeman Doppler maps are plotted in the right column. Both sets of potential 
field maps and Zeeman Doppler maps fit the model spectral dataset to reduced = 1-0 



3.1.1 



Bad phase coverage 



(1999), It is the p redicted potential fiel d for AB Dor in 1996 



In the case of high inclination stars, bad phase coverage leads 
to an increased amount of ambiguity in ZD maps, especi ally 
betwee n low-latitud e radial and meridional field vectors ( |Do-| 
nati & Brown 1997). We test the effects of bad phase cov- 
erage on potential field reconstructions and compare them 
with ZD maps for the same image and line parameters as 
those listed above but for a dataset with only 50% phase 
coverage. The potential field and ZD maps for this poorly 
sampled dataset are shown in Fig|5| As before, all recon- 
structions fit the observed dataset to a reduced = 1-0. 

We would expect the amount of cross-talk between field 
vectors to be much reduced in the potential field maps com- 
pared with ZD maps. The reason for this is that the line-of- 
sight component is strongest approximately 45° away from 
the centre of the stellar disk Hence, when the radial field 
contribution is reduced 40°-50° away from the disk centre 
due to poor phase coverage, the strong azimuthal field vec- 
tor contribution serves as a more accurate constraint in the 
potential field reconstructions. As the maps show in Fig. |[ 
we find that this is indeed the case. The potential field maps 
reconstruct a reduced amount of flux but retain the correct 
polarity and show less cross-talk than the ZD maps. This is 
also reflected in the calculations shown in Table |l|. 

3.2 Complex images 

We present potential field reconstructions to a more com- 
plex, realistic field distribution in this sectio n. The input 
image (Fig. ^) is of a potential field taken from Jardine et al 



December 23-25. Jardine et al. (1999) produced this image 
using the ZD radial field map from that epoch and a code 
originally develo ped to study the formation of filament chan - 
nels on the Sun (van Ballegooijen, Cartledge & Priest 1998). 
The model spectra generated using these maps are plotted 
along with the potential field reconstruction spectral fits in 
Fig. I 

Fig. ^ shows that the input and reconstructed poten- 
tial field maps below 75° latitude show good agreement. At 
9 > 75°, the discrepancy between input and reconstructed 
potential field maps increases. ZD maps also reconstruct 
more fiux near the poles than is present in the input maps 
but the disrepancy is not as great as with the potential field 
map reconstructions. At lower latitudes, on the other hand, 
it is clear that the potential field maps recover more accu- 
rate meridional field information compared to conventional 
ZD maps. 

It is worth noting that the polar region over which po- 
tential field maps recover the greatest spurious structure 
only represents a small area on the stellar surface. The ad- 
dition of this spurious structure is most likely due to the 
definition of the weighting scheme (Eqn. However, the 
accuracy of the potential field maps in the mid to low- 
latitude regions clearly demonstrates that the potential field 
code recovers more information about the surface field than 
is possible using conventional ZDI. The images produced 
by different weighting schemes (Fig. ^ show that different 
weighting schemes redistribute the recovered flux into differ- 
ent bin sizes on the surface. The overall field pattern is not 
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PF radial field map ZD radial field mop 




longitude longitude 

Figure 5. Maps showing the effects of poor phase coverage. Potential radial, azimuthal and meridional field maps are plotted in the 
left column and ZD maps are plotted in the right column. As before, all maps fit the spectral dataset to a reduced = 1.0. Observation 
phases are represented using tick marks. 




Velocity km/s 



Figure 4. Model Stokes V spectra generated using the input im- 
age shown in Fig. 3 are represented by the error bars. Spectra 
are plotted in order of increasing phase in steps of 0.05 phase. 
The predicted profiles from the potential field reconstructions are 
shown by the solid lines. The Zeeman Doppler predicted line pro- 
files are not plotted here but fit the model spectra to a similar 
accuracy. 



Table 2. Reduced chi-squared values evaluating the agreement 
for the reconstructions to the complicated input images. The first 
column lists the reconstructed set of images used for the calcula- 
tions; PF-1 stands for weights where n = 1; PF-2 means recon- 
structions where weights n = 2. The rest of the columns indicate 
the level of agreement between reconstructed images for each field 
orientation. 



Rec. ima 


ye 


■^rad 




Xme?' 


Complex 


ZD 


2.18 


2.21 


1.61 


Complex 


PF-2 


2.83 


1.94 


1.18 


Complex 


PF-1 


3.17 


2.25 


1.52 



changed but does affect the level of fit to the input PF map. 
The reduced measurements for these PF reconstructions 
(Table |) are clearly affected by the spurious structure at the 
pole in the radial field maps. However, they do show that in 
the case of reconstructions using the area density weighting 
scheme (n = 2) the azimuthal and meridional field recon- 
structions are more accurate. 



4 DISCUSSION 

Many limitations in ZDI arise from the sensitivity of cir- 
cularly polarized profiles (Stokes V parameter) to the line- 
of-sight components of magnetic fields. However, until spec- 
tropolarimetric instrumentation is sufficiently advanced so 
that all four Stokes parameters can be observed simulta- 
neously for cool stars this problem is unlikely to be re- 
solved. Currently, the four Stokes parameters can only be 
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Input rodial field map ZD radiol field map 




S 30 63 90 120 150 180 210 240 270 300 330 360 ' ' 33 60 93 120 150 180 210 240 270 303 330 363 

longitude longitude 



Figure 6. Input potential surface radial, azimuthal and meridional field maps calculated using Zeeman Doppler radial field map for AB 
Dor in 1996 December (Jardine et al. 1999). Potential field reconstructions are in the middle column and Zeeman Doppler maps are in 
the right column. Greyscales represent the flux in Gauss and the tick marks denote the observation phases used to generate model Stokes 
V spectra as before. 



measur ed in bright, highly mag netic A and B stars ( Wade 



et al. 2i)00a 



Wade et al. 2000b) 



The lack of physical realism in ZD maps also allows 
for unrealistic magnetic field distributions. By limiting ZD 
reconstructions to potential field distributions we can over- 
come some of these problems. In mapping potential fields, 
we introduce a mutual dependence on all three surface field 
vectors thereby recovering more field information. The ben- 
efits of this method are threefold; firstly, in regions where 
one vector is suppressed due to the inclination of the star 
(such as regions of low-latitude meridional field) there is still 
a reasonably high contribution from the other two and thus 
field information is reconstructed more reliably; secondly, in 
areas of bad phase coverage, azimuthal fields can sometimes 
still be recovered (as their signatures are strongest about 
45° away from the previous point of observation) and the 



mutual dependence of the field vectors can be exploited to 
recover some radial and meridional field information; and fi- 
nally, in regions such as starspots where the field is thought 
to be mainly radial but where there is little flux contribu- 
tion from all three vectors, the field information in the sur- 
rounding spot areas can be used to recover more of the lost 
polarized signal. 

In addition to these improvements, potential field con- 
figurations should allow us to test if the unidirectional band 
of azimuthal field seen in ZD images of AB Dor is necessary 
to fit the data. Such a band is incompatible with a potential 
field and its presence challenges our current understanding 
of magnetic field generation. 

If datasets are available that span a sufficiently long 
time-span it will be possible to monitor which modes are 
consistently active on these stars which will give us insight 
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PF radial field map - wl PF radial field map - w2 




30 60 90 120 150 180 210 240 270 300 330 360 ' ' 30 60 90 120 150 180 210 240 270 300 330 360 

longitude longitude 



Figure 7. Reconstructed potential radial, azimuthal and meridional field maps using diflferent weighting schemes. The first column shows 
reconstructions using weighting scheme with n = 1 and the second one shows maps reconstructed using n = 2. Greyscales represent the 
flux in Gauss and the tick marks denote the observation phases used to generate model Stokes V spectra as before. 



into the dynamo mechanisms operating in cool stars. Ex- 
trapolating these surface potential field models to the corona 
will enable us to model the coronal topology of the star and 
thus to predict where stable prominences may form in these 
stars. 



Finally, special care must bo taken when applying this 
technique to rotationally broadened Stokes V spectra from 
cool dwarfs. The weighting scheme used will have to be 
tested more thoroughly in order to ensure that we are not 
sacrificing realism for the "simplest images" (in terms of in- 
formation content). More realistic images may be obtained 
by incorporating a form of entropy which pushes the field 
distribution to a power law form, a possibility we are inves- 
tigating further. 
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